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Abstract 

We study a restricted solid-on-solid (RSOS) model involving deposition and evaporation with 
probabilities p and 1 — p, respectively, in one-dimensional substrates. It presents a crossover from 
Edwards- Wilkinson (EW) to Kardar-Parisi-Zhang (KPZ) scaling for p ~ 0.5. The associated KPZ 
equation is analytically derived, exhibiting a coefficient A of the nonlinear term proportional to 
q = p — 1/2, which is confirmed numerically by calculation of tilt-dependent growth velocities 
for several values of p. This linear A — q relation contrasts to the apparently universal parabolic 
law obtained in competitive models mixing EW and KPZ components. The regions where the 
interface roughness shows pure EW and KPZ scaling are identified for 0.55 < p < 0.8, which 
provides numerical estimates of the crossover times t c . They scale as t c ~ A~^ with <j) = 4.1 ±0.1, 
which is in excellent agreement with the theoretically predicted universal value = 4 and improves 
previous numerical estimates, which suggested <f> « 3. 
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I. INTRODUCTION 



The competition between different growth mechanisms is a characteristic of many real 
processes and has been the subject of intensive investigation in the last years. Many authors 
considered competitive growth models in which different dynamic rules are randomly chosen 
for the aggregation of the incident particles and applications to 

real systems were suggested . Il2j|. Such simplified models may mimic, for instance, the 
effects of large energy distributions of the incident atoms, which lead to different dynamic 
behavior as they arrive at the film surface. They usually show crossover effects from one 
dynamics at small times t or short length scales L to another dynamics at long t or large L. 

In many cases, a crossover from the Edwards- Wilkinson (EW) dynamics to Kardar- 
Parisi-Zhang (KPZ) growth Q| is observed. The Langevin-type equation 



dh 
dt 



v^ 2 h+ ^(V/i) 2 + v(x,t), 



(1) 



known as KPZ equation, is a hydrodynamic description of kinetic surface roughening, where 
h is the height at the position x in a <i-dimensional substrate at time t, ^ rep resents a surface 
tension, A represents the excess velocity and rj is a Gaussian noise ll4 MJ with zero mean 
and variance (rj (x, t) rj (x', tj) = D5 d (x — x'j 5 (t — t'). When A = in Eq. (JTj), we obtain 
the (linear) EW equation. Thus, if A is very small, the features of EW growth are expected 
at small times, and a crossover to KPZ behavior is observed at a characteristic time t c , when 
the macroscopic properties are affected by the overall nonlinear character of the process. In 
this paper, we will analyze universal and nonuniversal features of this crossover in lattice 
models through analytical and numerical methods. 

The roughness (or interface width) W(L,t) is the simplest quantity that indicates 
crossover effects. In lattice models, it is defined as 

1/2 v 



W(L,t) 



(2) 



for deposition in a (^-dimensional substrate of length L {hi is the height of column i at time 
t, the bar in h denotes a spatial average and the angular brackets denote a configurational 
average). In a typical EW or KPZ system, it scales for small times as 



W ~t 



(3) 
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However, when the crossover EW-KPZ is present, the roughness exhibits two growth re- 
gions, characteristic of EW and KPZ scaling (Pew < Pkpz in any dimension), as shown 
qualitatively in Fig. 1. At long times, the roughness saturates as 



W sat ~ L a . (4) 

t x is the crossover time to the steady state or saturation regime, also shown in Fig. 1. 

From plausible scaling arguments (reviewed in Sec. Ill), several authors suggested that, 
in d — 1, the EW-KPZ crossover takes place for small A at 

t c ~ A'*, (5) 

with a universal crossover exponent = However, to our knowledge the 

best known numerical estimate of this crossover exponent is ~ 3. It was obtained by 

n n 

Guo, Grossman and Grant |17( and by Forrest and Toral [19] through numerical solutions 
of the KPZ equation and data collapse methods. Recent works on lattice models confirmed 
the expected scaling relations for the growth and saturation regimes of KPZ, even in the 
presence of the EW-KPZ crossover p, but they were not able to improve the results for the 
EW regime (t <C t c ) or the crossover regions (t ~ t c ). Thus, neither a numerical confirmation 
nor a thoroughly justified refutation of the universality of the exponent = 4 was presented 
yet. 

On the other hand, an universal relation between the coefficient A and parameters of 
competitive lattice models with the EW-KPZ crossover was recently proposed by Braunstein 
and co-workers [?], Isj]. They considered processes where the aggregation of incident particles 
followed the rules of a KPZ lattice model with probability p and the rules of an EW model 
with probability 1 — p. The most studied representative jj, |f| is the competitive model 
involving ballistic deposition (BD - KPZ class) jstj and the Fami 



y model, also known as 



random deposition with surface relaxation (RDSR - EW class) 21]. The derivation of the 
corresponding KPZ equation from the stochastic rules of this class of models gives A ~ p 2 
for small p and is confirmed numerically for the RDSR-BD model p . 

In the present paper, we will study analytically and numerically a lattice model with 
the crossover EW-KPZ in d = 1, which is helpful to clarify the universal and nonuniversal 
relations in this crossover. The model is a restricted solid-on-solid (RSOS) one 



22 1 , in which 



deposition and evaporation of particles compete with probatilities p and 1 — p, respectively. 



EW behavior is found for p = 1/2, and KPZ behavior for p 7^ 1/2. We will derive analitically 
the KPZ equation for this process, which exhibits A~g = p— 1/2, where q represents a 
small relative probability of KPZ growth in the crossover region (p pa 1/2, q pa 0). This 
linear relation between q and A is confirmed numerically, and contrasts to the parabolic 
law found in other competitive models. Consequently, the A — p relation in the EW-KPZ 
crossover is clearly a model-dependent feature and not an universal law. On the other hand, 
our numerical work will also provide an estimate of the crossover exponent which agrees 
with the theoretically predicted universal value = 4, improving previous estimates which 
failed to confirm that prediction. This exponent is obtained from the scaling of t c (q), which 
is estimated from the intersection of the EW and KPZ behaviors, as illustrated in Fig. 1. 
The inherent difficulties of the numerical work, combined with the relatively simple, linear 
A — p relation, explain why estimating the crossover exponent is usually so hard. 

The rest of this work is organized as follows. In Sec. II, we will define precisely the 
discrete model, analitically derive its associated KPZ equation and confirm numerically the 
A ~ q relation. In Sec. Ill we will review the scaling arguments predicting = 4 and show 
the details of the numerical analysis which gives = 4.1 ± 0.1. In Sec. IV we summarize 
our results and present our conclusions. 



II. THE DISCRETE MODEL AND THE ASSOCIATED KPZ EQUATION 

In our competitive model, the deposit obeys the RSOS condition at any time, i.e. the 
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maximum height difference between neighboring columns is equal to the particle size a 
In the simulations, we consider a = 1. At each step of the process, a column of the deposit 
is randomly chosen. Subsequently, deposition and evaporation attempts are chosen with 
probabilities p and 1 — p, respectively. When evaporation is chosen, the top particle of 
that column is removed if the RSOS condition is satisfied after evaporation, otherwise this 
attempt is rejected. When deposition is chosen, a new particle is deposited at the top of that 
column if the RSOS condition is satisfied, otherwise this attempt is rejected. The time unit 
r corresponds to L attempts of evaporation or deposition in a substrate with L columns. In 
the simulations, we considered r = 1. 

This model was previously studied numerically by Amar and Family j^J, in order to 
show the universality of scaling functions and amplitude ratios for KPZ processes in d — 1. 
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However, that analysis was restricted to p > 0.75, consequently far from the region of 
EW-KPZ crossover. 

Now we construct the associated KPZ equation of this process starting from the master 



24j . following the standard method 



equation and performing a Kramers-Moyal expansion 
used in Refs. QBE I- 

First consider the deposition process according to the RSOS condition. The transition 
rate W(H, H') from the height configuration H = {hi} to the configuration H' = {h^} for 
this process is 

W(H, #') = -£ w^5 (h' k , h k + a) J] 6 fa, h 3 ) , (6) 

k jj^k 

where the 5 functions represent the condition that only the height of the column of incidence 
can be increased and w k describes the condition for aggregation: 

«4 0) = e (h k+ i - h k ) e (h k -! - h k ) , (7) 

where the Q(x) is the unit step function, defined as 0(x) = 1 for x > and Q(x) = for 
x < 0. Consequently, the first and the second transition moments are 



k\ 1] = £ (K - hi) W (H, H') = -0 (h i+1 - h) 6 (h^ - K) (8) 

T 



j {h- i -h t )W(H,H') = ° 

H 

and 

a 2 



K\- ] = £ (K - h) -hj)W (H, H') = -0 (h i+1 - hi) (hi—! - hi) 5(i,j). (9) 



'(2) _ 

H' 

For RSOS evaporation, the transition rate and the transition moments are those of Eqs. 
((ZJ), (jHJ) and © with opposite signs in the arguments of the functions. 

The Kramers- Moyal expansion of the master equation for the process provides the 
stochastic equation [24] 

dhi „ (1) ^_/7^ 



i-^^ + EV^V (10) 

where rjj is a Gaussian white noise with zero mean and co- variance (Vi(t)Vj{^)) = j)${t~ 
t'). For the competitive model, we obtain 

-± = p-Q (h i+1 - hi) (hi-! - hi) - (1 - p)-Q (hi - h i+1 ) (hi - hi- t ) + Di V i, (11) 

at r t 
where is constant. 
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In order to pass from the discrete description of the model to its continuum limit, we can 
use some analytical representation of the step function, which works in some limits. Many 
regularization for the theta step function have already been suggested, such as the hyperbolic 
tangent function 2^1 , and maximum function 2(| . This representation is expanded in Taylor 
series, so that 

9(x) = Co + ci x + c 2 x 2 + . . . (12) 

Inserting this expansion in equation (fTTj). we get 

dh; a r 



dt r 



c + ci(h i+1 - hi) + c 2 (h i+1 - hi) c + Ci(/i;_i - hi) + e 2 (/i;_i - K) 



- (1 -p)- c + cx(hi - h i+ x) + c 2 {hi - h i+ if c + ci(/ij - hi-i) + c 2 {hi - hi- X ) 
t 1 J L 

+ D iVi (13) 

In the continuum limit, a — > 0. In this limit, ac\ tends to a finite, nonzero value, since 
the angular coeficient C\ in the expansion of the theta function (Eq. IT^j) is of order 1/a. 
Moreover, a/r —>■ const, since this is the random growth velocity. We replace hi(t) by a 
smooth function h(x,t), whose coarse-grained derivatives are 



h i+1 - 2hi + hi_ x ~ a 2 V 2 h(x) 
h i+ i — hi ~ aVh(x) 



(14) 



Substitution in Eq. (fTB^) gives 



dh 



(2p - l)-c 2 + — c CiV 2 /i(a;) + 



+ (1 -2p)— (cl-2c c 2 ) \ Vh{x)\ 2 + 0{a 5 )+r]{x,t) 



(15) 



This equation must reproduce correctly the random deposition model, when all interac- 
tions between columns are turned off. Consequently, the best choice is to put c = 1. Also, 
all the above mentioned choices for representing the theta function, such as the hyperbolic 
tangent, lead to c 2 = 0. This is typical of odd functions, such as f(x) = 9(x) — 1/2. Thus 
we get 

^(f, t) = (2p-l)- + Cl -V 2 h + (1 - 2p) Cl 2 -(Vh) 2 + V (x, t), (16) 

Ot T T T 

which is the KPZ equation associated with the RSOS model with deposition and evaporation. 
All terms in the right hand side of Eq. (JTSJ) are finite quantities because a/r, ac±, Vh and 
aV 2 h are expected to have the same order of magnitude. 



It is interesting to recall that the choice of the value of 9(0) is arbitrary because the step 
function is nonanalytic at the origin. Thus, if our choice were 8(0) = 1/2, instead Eq. ((7j) 
we would have to represent the aggregation condition as 

= Q(hk+i-h k )Q(h h -i-h k ) [1 + 8(h k , h k + 1) +5(h k , h k -l) + 5(h k , h k + l)5(h k} h k - 1)] 

(17) 

where S(i,j) is the discrete Kronecker delta. As expected, this also gives the KPZ equation 
for the model, but the choice #(0) = 1 is suitable to represent the aggregation rule in a 
concise form. 

Comparison of Eqs. (|TB|) and shows that A varies linearly with q = p — 1/2. As 
expected, A < when deposition is dominant, and A > for dominant evaporation. Such 
inear relation is similar to that predicted for a single step model by Derrida and Mallick 
28| through a mapping into a one-dimensional asymmetric exclusion model. On the other 
hand, it contrasts to the A ~ p 2 law obtained by Muraca et al 171 for pure KPZ models (with 

n 

finite A) competing with pure EW models, such as the RDSR-BD model [6]. 

This linear A — p relation was confirmed numerically. The coefficient A can be calculated 
from the tilt-dependent growth velocity in the KPZ regime. If a given KPZ process takes 
place on an infinitely large substrate of inclination u, then A is related to the growth velocity 



/• as [29, 



30 



31] 



A = M (18) 



d 2 v 

' u=0 

(this form applies to d = 1, but is straightforwardly extended to higher dimensions). Several 
probabilities in the range 0.55 < p < 1 were considered for the simulations in substrates of 
length L = 10 4 . For each p, inclinations from u = 0.1 to u — 0.8 were considered, and the 
deposit was grown until times sufficient long for the KPZ regime to be attained. Average 
values were taken over 100 realizations for each p and u. 

Fig. 2a illustrates the method to calculate A from the growth velocities for three different 
values of p. The parabolic fits accurately represent the data behavior for all inclinations. 
Using those fits and Eq. (fT%|) . we obtained estimates of A for each p. In order to check 
the accuracy of these estimates, we also calculated the ratio (v — Vq) / u 2 (t> is the growth 
velocity at zero slope), and extrapolated that ratio to the limit u — > 0. The estimates of A 
agreed with those obtained from the parabolic fits within error bars. In Fig. 2b we show 
A versus q, which confirms the linear relation between those quantities for a large range of 
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values of g, in agreement with the KPZ equation obtained for the process. 



III. NUMERICAL STUDY OF THE EW-KPZ CROSSOVER 

First we recall the arguments that lead to the prediction of a crossover exponent = 4. 
In the works of Grossmmann, Guo and Grant (GGG) Il7f and of Nattermann and Tang 

n 

(NT) 18] ( see also [19]), this result is derived from multiscaling relations for systems with 
crossover EW-KPZ in d — 1. They proposed relations in the form 

W(L,t) = L°f(±-,£j, (19) 

where £ c ~ t l J z Ew j s a crossover length and zew — 2 is the dynamical exponent of EW 
processes. Assuming that t c scales as Eq. (JSJ), they obtained = zew/(&ew + z ew — 2) 
using scaling arguments. In d — 1, we have Pew = 1/4 and ccew = 1/2, which gives = 4 
in d = 1. This was confirmed by one-loop renormalization group calculations by NT. 

The same result also follows from the expected roughness scaling of KPZ in d = 1. 
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Assuming dynamic scaling in the nonlinear and saturation regimes, Amar and Family 
showed that the roughness scales as 

W(L,t) ^ L^g (\X\^j , (20) 

where g is a scaling function and where the dependence of W on the parameters v and D 
of Eq. Q was omitted. In the growth regime, it gives W ~ | A 1 1//3 t !/ 3 ((3 KPZ = 1/3 in Eq. 
|HJ). Now assuming that the crossover EW-KPZ takes place when the EW roughness (Eq. El 
with (3ew = 1/4) matches that of KPZ, as illustrated in Fig. 1, we obtain X 1 ^ 3 t c 1 ^ 3 ~ i c 1//4 , 
from which = 4 also follows. 

This last argument is the basis for the numerical calculation of t c using the roughness in 
the EW and KPZ regimes. First, it is necessary to calculate scaling amplitudes not shown 
in Eq. (JHJ): for the EW regime we have 

W E ~ At lf \ (21) 

and for the KPZ regime we have 

W K w Bt 1/3 . (22) 



S 



Matching these forms at t c we obtain 



.4 
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f c «(-) • (23) 



Consequently, t c can be determined from the estimates of the amplitudes A and B. 

Simulations of the RSOS model with deposition and evaporation were done for several 
values of p, in lattices with L = 10 5 , up to times approximately 10 6 . 100 deposits were 
generated for each p. 

In Fig. 3 we show W/t 1 ^ for small times t with p = 0.7 and p = 0.6. That ratio is 
expected to be constant in the EW regime. A narrow region 20 < t < 40 with this feature 
is observed for p = 0.7, while a wider EW region is found for smaller p. Here it is important 
to recall that other competitive models fail at this point because a clear EW region is found 
only for very small p, where the KPZ regime becomes difficult to be attained in simulation; 
one example is the model involving ballistic deposition and the Family model studied in Ref. 



The calculation of amplitude B is slightly more difficult because the ratios W/t 1 ^ 3 are 
not constant inside a time window long enough to extend to the maximum simulation times. 
In other words, the presence of significant corrections to scaling in Eq. ()22|) has to be taken 
into account. This can be done with the extrapolation of W/t 1 ^ 3 as a function of l/t 1//3 , 
as shown in Fig. 4 for p = 0.7 and p = 0.6 (see also Ref. jf|). Although the range of the 
variable 1/t 1 ^ 3 (abscissa of Fig. 4) is relatively small, it comprises almost two decades of 
the largest values of t. Good linear fits of the data are obtained in these large time regions, 
which suggests constant (but large) subdominant corrections to Eq.()22|). The amplitude B 
is estimated from the intersection of those fits with the vertical axis (t — > oo). 

For fixed p, different ranges of the variable 1/t 1 ^ 3 were chosen for the extrapolation of the 
data and the calculation of error bars in the estimates of the amplitude B. This procedure 
provides reliable and accurate final estimates of that amplitude. For instance, for p = 0.7 
(Fig. 4), we obtain B = 0.330 ± 0.004. We also observe that, while the amplitude A slowly 
varies with p (nearly 10% from p = 0.55 to p — 0.8), the amplitude B has a remarkable 
dependence on p. 

The estimates of t c obtained from Eq. (|2*3*|) are shown in Fig. 5 as a function of q = p— 1/2. 
Linear fits of different subsets of those data give 



The linear relation between q and A implies = 4.1 ± 0.1, which is in excellent agreement 
with the theoretically predicted value. 

Here it is important to recall that, in other competitive models such as the RDSR-BD one 

(|, the calculation of t c with accuracy was not possible. For instance, a clear EW growth 
regime (with Pew — 1/4) is observed in that model only for very small p, but in these 
conditions the KPZ growth regime (with Pkpz = 1/3) is not attained within a reasonable 
simulation time. This may be a consequence of the typically huge scaling corrections of BD 

3^ . However, we believe that the main reason for those difficulties is the parabolic A — p 
relation, which significantly reduces the range of p where both regimes can be numerically 
analyzed. 



IV. CONCLUSION 



We studied a competitive growth model in 1 + 1 dimensions involving RSOS deposition, 
with probability p, and RSOS evaporation, with probability 1 — p. This model may be 
viewed as a discrete realization of the continuum KPZ equation with an adjustable nonlinear 
coupling A related to p. Its corresponding KPZ equation is derived, showing that |A| linearly 
increases with q = p — 1/2, so that the process belongs to the EW class for p = 1/2. This 
result is confirmed numerically by calculation of tilt-dependent velocities for several values of 
p. It contrasts to the parabolic A— p relation obtained for competing models involving a KPZ 
and an EW process, which shows that this relation, although being of wide applicability, is 
not universal. 

We also calculated numerically the scaling amplitudes of the EW and KPZ growth regimes 
for several values of p. From these quantities, estimates of the crossover times t c were 
obtained. They scale as Eq. (j5J) with = 4.1 ± 0.1, in excellent agreement with the 
theoretically predicted value of the crossover exponent. This result improves previous ones, 
which suggested <p ^ 3 from simulations of the KPZ equation. We believe that this work 
provides an important, possibly definite confirmation of scaling relations predicted for the 
EW-KPZ crossover in d — 1. 
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FIG. 1: Typical time evolution of the roughness of a system with an EW-KPZ crossover at time 
t c and saturation time t x . 
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FIG. 2: (a) Growth velocity v as a function of inclination u for the competitive model with p = 0.55 
(triangles), p = 0.65 (crosses) and p = 0.75 (squares). Dashed lines are parabolic fits of each set of 
data, (b) Estimates of A as a function of the reduced probability q = p — 1/2 (squares) and a least 
squares fit of the data (dashed line). 
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FIG. 3: W/t 1 ^ at small times t for the competitive model with p = 0.7 (circles) and p = 
(squares) in a large lattice (L = 10 5 ). The inset shows a zoom of the data for p = 0.7 in the 
region. The dashed lines are linear fits of the data in those regions. 
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FIG. 4: W/t 1 / 3 versus 1/i 1 / 3 at long times, for p = 0.7 (circles) and p = 0.6 (squares) in a large 
lattice (L = 10 5 ). The dashed lines are linear fits of the data. 
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FIG. 5: Crossover time t c versus q = p — 1/2 for the competitive model with 0.55 < p < 0.8. The 
dashed line is a linear fit of the plot, with slope near —4.1. 
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